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INTRODUCTION 


Various forms of quartz have long been studied and constitute 
an important body of information for geologists and geophys- 
icists. Quartz is a ubiquitous mineral. It is a major con- 
stituent of many sedimentary rocks and an important component 
in some igneous and metamorphic rocks . The polymorphism, 
stability and reaction kinetics of quartz in several of its 
six forms is important in concepts of the mineralogy and physi- 
cal properties of the earth's interior. An example is some of 
the several possible stable phase assemblages or transformations 
of feldspar between 30 and 200 kb: 

1. Oligoclase jadeite + gossularite + 
kyanite + Si 02 (quartz, coesite, stis- 
hovite) . James (1969) . 

2. Mi crocline -»■ hollandite. Ringwood, 

Reid and Wads ley (1967) . 

The density of the first assemblage depends upon the thermo- 
dynamic state of quartz. In turn, these phase assemblages 
might be used in more elaborate models of the mantel which 
must be consistent with density variations within the mantle. 

This report describes and presents the results of a study of 
shock wave loading of Coconino sandstone from Meteor Crater, 
Arizona. These results extend and add to previous work des- 
cribed by Jones, et al. (1968); Ahrens and Gregson (1964); 
and Gregson and Grine (1963) . Combining the results of these 
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workers , we attempt to describe the polymorphism of quartz 
(i.e., a-quartz, coesite, stishovite) under high rates of 
loading at high pressures. The hugoniot of sandstone leads 
us to believe that mixed regions of both a-quartz/coesite and 
a-quartz/stishovite do exist at high rates of loading not only 
in sandstone but also in crystalline quartz as well. 

Experimental techniques employed in the present work to achieve 
high pressures are described in detail elsewhere (Jones, et al., 
1968 and Shipman, et al., 1968). The technique involves the 
use of an accelerated reservoir light gas gun to launch a flat 
plate of one of several metals with well-documented shock wave 
properties. Metal plates are made to strike an instrumented 
specimen and the data parameters of the shock wave are recorded. 
This technique allows the generation of pressure states higher 
than previously attained through the conventional use of explo- 
sives because impact velocities as high as 8 km/sec are obtained. 

Conversion of recorded data to pressure-volume states in the 
shock front is possible by using precise measurement of the im- 
pactor velocity and shock wave velocity in the specimen together 
with the known hugoniot equation of state of the impactor. 

The Rankine-Hugoniot conservation equations are then used to 
calculate the pressure-volume and energy state behind the shock 
wave front in the sample. 


SECTION I 


SHOCK WAVE DATA 


The physical properties of Coconino sandstone tested in this 
study are presented in Table I . The general petrologic des- 
cription and physical properties of this Coconino sandstone 
are nearly identical to those reported by Ahrens and Gregson 
(1964) . Measurements of the acoustic velocities using a 
pulse transmission overlap technique are included for further 
completeness of material characterization. 

A total of 17 experiments were performed using an accelerated 
reservoir light gas gun facility. Measurements of hugoniot 
shock states in Coconino sandstone ranged from about 40 to 
1400 kb. This range of shock states was achieved by varying 
the launch velocity of the flat plate impactor from approxi- 
mately 1 to 8 km/sec and using impactors of varying shock im- 
pedances. The experimental techniques will be discussed below 
in aspects related to data comparison. 

Data presented in Table II are also displayed in Figures 1 and 
2 showing shock velocity-particle velocity and pressure-specific 
volume relationships respectively. Data tolerances and error 
bars are included to illustrate the scatter recorded at each 
represented shock state. 

The simplest representation of hugoniot data as a description 

of behavior in shock velocity-particle velocity coordinates 

(U , u ). A least squares analysis is used to fit a smooth 
s p 
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TABLE I 


PHYSICAL PROPERTIES OF 
COCONINO SANDSTONE 


Mineral Cbmposition 

97% Quartz 
3% Feldspar 

Trace Clay and Heavy Minerals 
Grain Size 


Average and Modal Grain Size - 0.117 - 0.149 mm 
Range of Grain Diameters 0.062 - 0.71 mm 

Density 

3 

Bulk Density - 1.98 gm/cm , 

Grain Density - 2.6 7 gm/cm"^ 

Porosity (Calculated) 


25% 

Strength 

Unconfined Crushing Strength normal to bedding 
3.14 X 10^ dynes/cm^ dry, 3.64 x 10^ dynes/cm^ 
saturated with water (block size 2.2 X 2.2 X 5 cm; . 

Sonic Velocities 

Longitudinal Cj^ = 4.188 km/sec ± .08 (Vacuum Dry) 

Shear C^ = 2.68 km/sec ± .3 (Vacuum Dry) 

Description 

Specimen is weakly to moderately well cemented with 
silica, in the form of quartz overgrowths on the grains. 
Subparallel laminae 5.0 to 17.5 mm thick are separated 
by thin laminae 0.5 mm thick containing more than average 
amounts of silt and clay sized grains. A petrofabric study 
of C-axis orientations can be found in Gregson and Grine, 
(1963) . 
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TABLE II 

HUGONIOT DATA FOR CONCONINO SANDSTONE 

= 1.98 ± .03 gm/cc 


Irapactor Material 
and Velocity 

Remarks 

Shock 

Velocity 

Particle 

Velocity 

Shock 

Pressure 

Specific 

Volume 

Shot No. 

(km/sec) 

Thickness in mm 

(km/sec) 

(km/sec) 

(Mb) 

(cc/gm) 


Cu 

.851 i 

No Shim 

2.95±.12 

.74 

.0431.002 

-.3781.005 

S-241 

Cu 

1.211 

Wedge Exper, 
2. 5-12. 5mm 

3.20±.06 

1.039 

.0661.001 

.3411.003 

S-257 

Cu 

1.588 

5.84 

3.674±.001 

1.330 

.()968 

.3222 

S-153 

Cu 

1.965 

5.84 

4.10 

1.63 

,132 

.304 ! 

S-156 

Cu 

2.054 

5.84 

4.269±.002 

1.700 

.1437 

.3040 

S-247 

Cu 

2.650 

5.84 

4.49±.10 

2.18 

.1941.004 

.2601.005 

S-1127 

Cu 

3.155 , 

5.84 

4.84±.02 

2.57 

.2461.001 

.2371.002 

S-152 

Cu 

3.528 

5.84 

5.25±.02 

2.86 

.2971.001 

.2301.001 

S-239 

Cu 

3.898 

5.84 

5.66±.01 

3.12 

.3491.001 

.2271.001 

S-154 

Cu 

4.063 

5.84 

5.79+.04 

3.25 

.3731.002 

.2221.005 

C-1128 

'Cu 

5.540 

5.84 

7.57+.03 

4.30 

.6441.003 

.2181.002 

C-1130 

Cu 

5.731 

5.84 

7.79±.14 

4.43 

.6841.011 

.2181.010 

C-1179 

Cu 

6.642 ; 

5.84 

8.72+.06 

5.07 

.8751.006 

.2111.002 

S-160 

Cu 

6.652 ' 

* 

5.84mm 

8. 821.03 

5.07 

.8861.003 

.2151.003 

S-158 

Cu 

6.662 ’ 

• 

4 . 6 3mm 

9.081.01 

5.05 

.9081.002 

.2241.001 

S-159 

Cu 

7.895 

* ! 

3.94mm 

10.09 

5.94 

1.187 

.208 

S-155 

FS- 

11 7.871 

• 

3.94mm 

11.20 

1 

6.43 

1.426 

.215 

S-157 


* Shim Thickness = .082mm Brass 
FS-Fansteel-77 
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Figure 1 


Shock Velocity versus Particle Velocity 
for GM Data on Coconino Sandstone 










curve to the U - u data. Most commonly, a linear function 

■ ' ' S ' ' p ■ ■ i' .. .... 

of the form + S u^ is employed as an auxiliary equation 

of state which allows interpolation "^between data. Where dis- 
continuities are evident, linear line segments are here em- 
ployed for the same purpose but only over the range of particle 
velocities of the segment. In general, but not necessarily, 
departures from linearity are evidence of changes in structural 
behavior e.g. plastic failure, melting, polymorphism in shock 
loaded materials (McQueen, et al. , 1967) . 

Comparison of the shock velocity-particle velocity and pres- 
sure-specific volume data of Ahrens and Gregson (1964) and 
this work is made in Figures 3 and 4 showing our interpreta- 
tion of these data. 

Region I represents the elastic behavior of sandstone established 
by Ahrens and Gregson (1964) where the HEL (Hugoniot Elastic 
Limit) was measured to between 4.1 to 5.1 kb. 

The next two regions use the data from the HEL (Hugoniot Elastic 
Limit) as determined by Ahrens and Gregson (1964) and lower data 
points reported here. 

REGION II 

U = 1.751 + 1.495 
s 

0.3<u <0.8 km/sec 
P 

REGION IIL 

U = 1.380 + 1.652 u km/sec 
s p 

0.95<u <1.8 km/sec Std. Error of U = 0.112 km/sec 
P s 


u km/sec 

p 

Std. Error of U = 0.103 km/sec 
s 
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SHOCK VELOCITY ' km/sec 



PARTICLE VELOCITY - km/sec 


Figure 3 Shock Velocity versus Particle Velocity 
Data of Sandstone 
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PRESSURE - Mb 
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T 
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SANDSTONE 



Figure 4 Hugoniot for Sandstone Data 
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REGION IV 


U = 3.303 + 0.588 u km/sec 
s p 

1.8<u <2.55 km/sec Std. Error of U„ = 0.126 km/sec 
P s 

REGION VA 

U = 1,172 + 1.427 u km/sec 
s p 

2.55<u <3.35 km/sec Std, Error of U = 0.02 km/sec 
p s 

REGION VB 

U = 0.422 + 1.661 u km/sec 
s p 

Up>3.35 km/sec Std. Error of = 0.132 km/sec 

REGIONS VA AND VB COMBINED: 

U = 0.567 + 1.633 u km/sec 
s p 

Up>2.3 km/sec Std. Error of = 0.117 km/sec 

Data for Coconino sandstone of Ahrens and Gregson (1964 ) , 
Gregson and Grine (1963) are employed here because their 
material description matches ours almost identically. The 

only material difference is the initial density ranging from 

3 3 

1.96 to 2.03 gm/cm versus the 1.98 ± .03 gm/cm average 

value measured for the present work. 

Use of the Rankine-Hugoniot jiimp conditions in this work and 
in previous works assumes a steady state behavior. In order 
to determine whether steady-state conditions existed in the 
present experiments, several tests were performed which would 
exhibit non-steady waves by a variation of average shock 
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velocities with thickness. One such test employed a wedge of 
sandstone which was impacted on one face and the arrival of the 
shock at the tilted rear surface was ,re corded using a high speed 
rotating mirror camera. The impact conditions produced a shock 
in sandstone of 66 kb, which was nearly 'Steady over the target 
thicknesses of 2.5 to 12.5 mm. Shock wave breakout as witnessed 
by an unshimed mirror resulted in a velocity of 3.20' km/sec. 

The slightly ragged cut off showed only dispersion of .06 km/sec 
('x.2%) across the wedge with the maximum velocities occurring be- 
tween 5 to 9 mm. Five to 9 mm is the thickness at which most of 
our hugoniot measurements were performed. 

The above velocity dispersion corresponds to about a .05 mm (1/2 
the average grain diameter) unevenness of the shock wave front 
which is equivalent to the surface roughness of the target speci- 
men. 

Tests were performed with and without' rear surface shims in an 
effort to assess their effects on recorded wave arrival times . 
Because particle velocity (or free surface velocity) is riot 
measured directly by our pin technique, any difference in this 
parameter affected by. shims would .appear as data scatter. The 
scatter of ^data about the hugoniot as interpreted above is 
shown by the stancJard errors to be relatively small. 

Three .other tests were performed at much higher shock stresses 
with the more sensitive pin techniques to check for non-steady 
propagation. Each target was a different thickness but the im- 
pact conditions were held nearly identical. These experiments, 
as seen in Figure 1 and Table I, display no consistent behavior. 
Their '^2% scatter is larger than the experimental uncertainty 
and is not in the order of descending velocity with increasing 
thickness or vice versa. Our. conclusion, is that the specimen 
or shot reproducibility is probably no better than 2%. Within 
this bound, shocks appear steady over the measured data. 
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SECTION II 

PRESSURE-VOLUME PHASE DIAGIU^ FOR QUARTZ 


The interpretation of the hugoniot data from the previous 
section is employed in the following proposed description of 
polymorphic behavior of quartz under shock loading. Breaks 
or discontinuities in the sandstone and solid quartz hugoniots 
together with static pressure-volume data are used to estimate 
a P-V phase diagram for quartz. Primary concern is directed 
toward conditions for existence of a-quartz, coesite, stisho- 
vite phases along tl^e hugoniot. 

Inference of polymorphic behavior of materials under shock 
(high strain rate) loading is usually made on the basis of 
breaks or discontinuities in the hugoniot or shock velocity- 
particle velocity plot. Recovered material, shock unloading 
paths, and correlation with static P-V and P-T data are also 
used, when available, as guides in the interpretation of the 
high rate loading data. 

Shock behavior of idealized models of a material undergoing 
a polymorphic transition is discussed in detail below. Four 
such models are illustrated in Figures 5. -and 6 . These models 
are based on the known behavior of quartz but do not include 
all possible transitions, endothermic reactions, melting, phase 
lines such that dP/dV>0, or those that are centered in a high 
temperature phase. 

Figure 5 illustrates the shock loading of a material which has 
a polymorphic transition that occurs in time intervals less than 
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(Figure 5b) and greater than (Figure 5c) the rise time of the 
shock front. The heavy lines denote phase boundaries. The 
dashed line path ABCD, (Figure 5a) denotes the static or equi- 
librium states of the material undergoing a phase transition. 





b. c. 


Figure 5 (a) Idealized Models of Polymorphic Phase 

Transformation (b) Having Transition Times 
Less Than and (c) Greater Than the Risetime 
of the Shock Front 
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Path ABEF (Figure 5a) is typical of a material in which the 

transition time of the phase transfoarmation is less than the 

rise time of the shock front. The area difference beneath 

paths ABCD and ABEF represents the non-equilibrium or "waste 

heat" energy between the static and dynamic loading states and 

results in a temperature increase at the shock front. This 

type of phase transition is described in some detail by Rice, 

McQueen and Walsh (1958) . The - u^ plot of a material 

having a hugoniot represented by path ABEF is shown in Figure 

5b. This is typical of a material in which a single stepped 

shock front occurs from A to B and again near F. A two-stepped 

shock front is stable in the region near B to E. Breaks or 

discontinuities in a U - u or P-V plot for this idealized 

s p 

model would be an adequate indication of the onset and comple- 
tion of a polymorphic transition providing other effects such 
as yielding can be accounted for. 

Path AGF (Figure 5a) is one of several representations that is 
possible for a material in which the phase transformation occurs 
in time intervals greater than the risetime of the shock front, 
the difference between path ABCD and path AGF again represents 
the non-equilibrium energy conditions between the static and 
dynamic loading states. The - u^ plot is shown in Figure 5c. 
This is typical of a few materials having a large precursor decay 
in the double stepped wave-front from B to G. Single stepped 
shock fronts are stable from A to B and near F. This model il- 
lustrates one of several ways in which "time-dependent phase 
transitions" can exhibit a transition from one phase to another. 
Horie (1966) discusses such effects in some detail as does 
McQueen, et al. (1967). Path AGF in Figures 5a and 5c can be 
viewed in two ways: a pressure in excess of B is needed to 
drive the transition which then by stress relaxation dilates 
to B within the duration of the shock experiment, or the tran- 
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sition time is greater than the risetime of the shock front 
(but less than the duration of the experiment) resulting in 
phase A being Overdriven and then relaxing. This type of be- 
havior can be indicative of the existence of a polymorphic 
transition but its onset cannot be determined until the de- 
caying double shock has traveled far enough into the material 
to reach a steady state propagation. 



a. 




b. c. 


Figure 6 (a) Idealized Models of Polymorphic Phase 

Transformation of a (b) Solid and (c) Porous 
Material Having a Transition Time Greater 
Than the Risetime of the Shock Front 
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Figure 6 illustrates two additional shock loading polymorphisms . 
Path AEFG (Figure 6a) is somewhat similar to path ABEF for Fig- 
ure 5a except that a two stepped shock cannot occur and a single 
shock will be stable. Again, the difference between paths ABCD 
and AEFG represents the non-equilibrium condition between the 
static and dynamic loading states. Point E may represent an 
overdriven state that could relax to point B. We shall see from 
our experiments that the amount of overdriving and relaxation 
is difficult to determine and that B and E are "close" for pres- 
sure levels encountered here. Thus, a P-V or shock U - u 

s p 

path AEFG (Figure 6a and 6b) could be an adequate indicator of 
phase boundaries . ' 

Path HIJK (Figure 6a and 6c) is our model of a porous material 
undergoing a polymorphic phase transition. Path H to I is the 
regime where pore collapse occurs and where shock temperatures 
are expected to increase rapidly. Complete collapse of pores 
is generally assumed to have taken place before I but this need 
not always to be the case. Differences between paths ABCD and 
HIJK represents the non-equilibrium energy conditions as before. 
Point I is assumed "close" to point B because it may lie to 
either side of A' A' because of the competing effects of the 
higher temperatures available to help initiate and drive the 
transition and time dependency of the pore reduction. Again 
a single shock will be stable throughout the transition region 
I to J. Single shocks are found in this case because the large 
initial volume of the porous material allpws the rayleigh line 
representing the "jump condition", to penetrate, to every part 
of the hugoniot. . 

Figure 7 shows where models in Figures 5c, 6b and 6c have been 
applied to the hugoniots of porous and solid quartz and illust- 
rate the polymorphism of quartz. (The. solid quartz hugoniots 
was abstracted from Wackerle (1962) and Fowles (1967).) 
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SOLID QUARTZ 


POROUS QUARTZ 




a 


b 


Figure 7 Interpretation of Solid and Porous 
Quartz Hugoniots 

Data from shock unloading, shock recovery experiments and static 
(low strain rate) experiments are used to generate a phase dia- 
gram for quartz (Figure 8) . The phase lines shown are located 
within limitations imposed by experimental errors of both the 
solid and porous hugoniot data. These boundaries delimit the 
P-V states at which mixed regions are totally converted for 
single phase and indicate only approximately the conversion 
temperatures. They establish nothing about the condition of 
mixing . 

Figure 8 is a pressure-volume phase diagram for quartz and 
summarizes our interpretation of its overall behavior at high 
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Figure 8 Pressure-Volume Phase Diagram for Quartz 
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pressures. The coesite to s.tishovite pJiase line is based on: 
the 'highest pressure-volume state (mean, normal .‘^pressures be- 
tween 360 and 385 kb) at which a change in slppe; occurs in 
the solid and porous hugoniots 'and-'a reberif paper by; Ahrens 
et al. (1970) using our data- to- esiimate- the sti^hovite iso- 
therm and by McQueen, et al. (1967) . The a-quartz-coesite 
phase line is based on the specific volume of ppesite at 1 bar 
by Sclar, et al. (1962) and on the indicated change in slope 
of the sandstone hugoniot at ^^154 kb. 

The a-quartz-coesite mixed region may be evident in the solid 
shock wave data as the extremely large stress » relaxation 

(Wackerle, 19 62; Fowles;, 1967; and Ahrens .and Gregson, 19 64) 

■'3' 

that takes place at or near V = .375 cm /gm which is commonly 
labeled the Hugoniot Elastic Limit (HEL).-. Ahrens and Duvall 
(1966) have made the most detailed study of this stress relaxa- 
tion in quartz at its HEL but the mechanism to produce this 
decay is an open question. Certainly time dependent generation 
of even a small amount of coesite (thus placing the material in 
the mixed a-quartz coesite region) could result in a redistribu- 
tion of shear stresses which could initiate the relaxation. 

The strongest argument for this mixed region is the break in the 
sandstone hugoniot. 

The previous behavior of quartz postulated by McQueen, et al. 
(1963); Ahrens and Gregson (1964); and Jones, et al . 1968), was 
that the stress relaxation observed in quartz is the result of 
an unknown mechanism and that partial transformation to coesite 
does not occur above the HEL. They went on to ‘postulate a 
direct transformation of a-quartz to stishovite at higher pres- 
sures. Coesite was not considered to exist along the hugoniot 
in these studies because of the extremely long time and high 
temperature required for the reaction to occur as observed in 
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static experiments. However, static phase transformations in 
calcite which also require long times and high temperature were 
found to occur within the thickness of the shock front by Ahrens 
and Gregson (1964) and Jones and Froula (1969) . This has been 
noted for other phases in other materials. Breed and Venable 
(1968) find large time variations in shock induced phase trans- 
formations . The great shear deformation forces produced in 
shock fronts is usually assumed to be the driving mechanism. 
Dremin and Breusov (1968) have recently reviewed this aspect 
of dynamic processes. 

Coesite is usually found as traces and small amounts in the 
recovered thetomorphic glass from meteor impact craters and 
nuclear detonations (Deribas, et al. 1966; Chao, 1967). 

Bunch, et al. (1968) also report shock induced structural 
disorder in quartz above 150 kb. Structural disorder might 
be evidence for partial coesite conversion followed by partial 
reversion of the crystal lattice by unloading. Deribas, et al. 
(1966) have published the only report claiming to find coesite 
in laboratory scale shock experiments above 150 kb. Many U.S. 
investigators have looked for coesite in recovered material 
from laboratory shock experiments but have found none 
(Chao, 1970) . 
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SECTION III 


TEMPERATURE CALCULATIONS 


In the preceeding section, static equilibrium pressure-volume 
states are compared with high-strain rate nonequilibrium pres- 
sure-volume states. A similar comparison can be made with 
pressure-temperature states. Pressure-temperature states for 
low strain rate equilibrium conditions are usually measured. 
Because shocks are of short duration ('^^)Jsecs) , temperature 
measurements are beyond the present state of the art. 

The low-strain rate P-T phase diagram for quartz is illustrated 
in Figure 9. References for the numbered boundaries are listed 
in Table III. The quartz-coesite boundary is obtained from 
Holm, et al. (1967); Bell, et al . (1965); Boyd, et al. (1966); 

Boyd and England (1960); Yasukawa (1963); Takahaski (1963); 
Kitahara and Kennedy (1964); Dac'hille and Roy (1959); MacDonald 
(1956) ; and Griggs and Kennedy (1956) . The a-quartz - 3-quartz 
boundary is obtained from Gibson (1928) , Yoder (1950) , and 
Cohen and Klement (1967) . The tridimite and cristobalite sta- 
bility fields are taken from Holm et al. (1967) ; Stishov (1963) ; 
Sclar, et al. (1962); Bell and Boyd (1968); and Wentorf (1962). 
In order to generate an equilivalent high strain rate P-T phase 
diagram for quartz, shock temperatures must be calculated. 

Calculation of temperatures in solid materials is usually done 
by methods first outlined by Walsh and Christian (1955) . 

The general problem of deriving a complete P-V-T equation of 
state from shock -data are considered by Cowperthwaite 
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(1966, 1969) and by Cowperthwaite and Ahrens (1967). Ahrens 
and Gregson (1964) , have made the first calculations in sand- 
stone and the present analysis follows their method. 



DEGREES KELVIN 


Figure 9 Low-Strain Rate Phase Diagram for Quartz 




TABLE III 


LOW STRAIN RATE DATA FOR EQUILIBRIUM PHASE 
DIAGRAM OF QUARTZ (FIGURE 9) 


1. Yoder (1950) 

Gibson (1928) 

Cohen and Element (1967) 

2. Tuttle and Bowen (1958) 

3. MacDonald (1956) 

Griggs and Kennedy (1956) 
Dachille and Roy (1959) 


4. Boyd and England (1960) 
Bell, et al. (1965) 

Holm, et al. (1967) 

Boyd, et al. (1966) 

Yasukawa (1963) 

Takahaski (1963) 

Kitahara and Kennedy (1964) 


5. Holm, et al. (1967) 
Stishov (1963) 

Sclar, et al. (1962) 
Bell and Boyd (1968) 
Wentorf (1962) 


6. Mosesman and Pitzer (1941) 


The agreement between the present calculation and that of 
Ahrens and Gregson (1964) is satisfactory. Error estimate 
based on possible variations of thermodynamic parameters 
made by Ahrens and Gregson (1964) indicate that these temper- 
atures are conservative. The same arguments also hold for 
the present calculations. We also assume that the Gruneisen 
ratio, to be volume dependent and tempertature independent, 
assuming only longitudinal lattice modes and neglecting shear 



lattice and optic lattice modes. These assumptions may re- " 
suit in incorrectly calculated temperatures (Anderson, 1968) 
but no data exists - to include these effects or to estimate-the 
errors ... -‘i'' ' ^ 


To Calculate temperature in sandstone, the hiigoniot internal 
energy per unit iriass is partitioned into compressional energy, 
Ej^ , • and thermal ehergi? , ; - ^ 


1/2 (P^ + P^)(V^ . = ®lc 


Where the o subscript is the initial state and the 1 subscript 
iS' the final state. V is the specific . volume state-and P is 
considered as hydrostatic pressure because of the initial small 
shear . stress , early crushup and assumed stress equalization 
under shock loading. The total; internal energy is obtained 
from the experimentally determined hugoniot; the compressional 
energy is determined from an isotherm. The area difference be- 
tween these curves is considered to be the thermal energy. 


- f'o' <''o - ''i' - =k 

o 

Where is the specific heat capacity at constant volume. 

The specific heat capacity at constant volume can be obtained 

from experimental measurements of C , the specific heat capa- 

P 

city at constant pressure, and from other thermodynamic para- 
meters such as a, the volume coefficient of thermal expansion 
and the isothermal bulk modulus through the Gruneisens Law. 
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C - C = a^TV^ (~~) 

p ^av^T 


Because (ap/av)^ is estimated for small increments along the 
hugoniot and below the Debye temperature is a function of 
temperature for each phase of quartz, temperature calculations 
are broken up into a series of regions where each of the above 
parameters can be handled with simplicity. Figure 10 shows 
the six regions which were used for temperature calculations 
in the shock loading of Coconino sandstone. Figures 11 and 12 
show the experimental specific heat capacity data as well as 
the empirical fits used in the calculations . 


For the elastic regime. Region I, the temperature is estimated 
by assuming adiabatic compression. This assumption is justi- 
fied because the elastic velocity is approximately equal to 
the adiabatic compressional wave velocity. Thus the hugoniot 
and the adiabat are close. The temperature increase is thus 
calculated as : 

r (V -V) 

T = T exp ^ 

o ^ V 

Where r is the Gruneisen coefficient and is assumed to be a 
constant . 


For sandstone in Region I, it is difficult to obtain consistent 
values of T. From specific heat at constant pressure and volume, 
r is (Kelley 1949) : 


r = (VCp) -1 

2 

V7here C = a + bT + cT , a, b , and c are empirical constants 
P 

and T is absolute temperature. 
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Figure 10 


Coconino Sandstone Hugoniot Showing Pressure 
Regions Used for Temperature Calculations 


TEMPERATURE- 


500 h- Kelley (I960) 


a to e Quartz 
Transition 


McSkimin.et al 
(1965: 


a Quartz Data Taken from 
Mosesman and Pitzer, (1941) 
McSkimin.et al, (1965) 

Janaf Tables (1965) 


Cp - Cal/Mole'"K 


Specific Heat Capacity as a Function of 
Temperature for a-Quartz, 6-Quartz, 
Tridimite and Cristobalite 







Using the definition of F as: 


r = aK/C p 

P o 

volumetric thermal expansion coefficient 
adiabatic bulk modulus 
specific heat at constant pressure 
initial specific density 

and using handbook values (Clark, 1966) , F has a value of 
0.01 or less. Using the smaller F produces an increase in 
temperature upon compression to the HEL (4.6 kb) of a few 
tenths of degrees Kelvin. 

Region II in sandstone is the regime between the HEL and the 
first break in the hugoniot. It is a region where the pores 
in the rock are being closed (points 0 to 1 in Figure 8) . The 
compressional energy, Ej^, is neglected because of a lack of 
experimental data in this pressure regime. The increase in 
surface energy of quartz grains during crushup has been esti- 
mated by Brace and Walsh (1962) but this energy is small com- 
pared to the compressional energy or shock energy and is also 
neglected for this region. Thus, all the shock internal energy 
is assumed to be converted into thermal energy or, 

rT, 

1/2 (P^ + P^) (V^ - V^) =l " C^(T)dT 

'^o 

An initial temperature, T^ , or 300° K is assumed in all calcu- 
lations of the shock temperature, T (neglecting temperature 

hi 

increase at the H.E.L.). 


Where 


a 

K 

C 



I 


Specific heat capacity at constant volume for a-quartz in this 
temperature regime is computed incrementally form: 


= Cp - T aVoP/3V)^. 

3 —1 

Where T(° Kelvin) and V (cm gm ) are average values of each in- 

-5 -1 

crement (9P/9V) along the hugoniot. A value of IX 10 °K 
is used for the volume coefficient of thermal expansion, a. 

Specific heat at constant pressure is taken from Figure 11 as: 


C = -.13244 + 0.0448T - 3.074 X 10 
P 

Where C has units of calories/mole °K and T is in units of 
P 

° Kelvin. C is adjusted for an assumed linear increase in 
P 

mass/unit volume from the 25% porosity condition of the sand- 
stone to that of solid quartz over this region. 

In Region III, a new energy partition was assumed because the 
porous material is considered to have been compressed to a solid. 
The internal energy was assumed to consist of a compressible and 
a thermal part. 


(P,-P^) (V„ - V, ) 


III 


i""o' ^''o " ''i' _ f 

2 J 


C (T)dT + E, 

V k 


Where is the energy of isothermal compression. 


Thus E, = 
k 


( 


V 


- PdV. 
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If the Murnaghan form of the isotherm is assumed: 


P = 



Where 3 = 42.75 kb and ? = 8.390 (Ahrens and Gregson, 1964). 


The equations are solved for and iterated on and 

until the temperatures agreed within a few tenths of a percent. 

Calculations were performed along the Murnaghan until 

3 

V = .355 cm /gm (point 2) was reached. This is interpreted 
as the mixed phase boundary between a-quartz and coesite regions . 
At this point, the energy from this isothermal compression using 
the Murnaghan form was held constant and energy was added in- 
crementally along the flat portion of the 300°K isotherm. 


V = 

.355 

.355 

- 3 

0 

\ 

!-• 

1 

dV + 

V = 

.375 J 

. 395 

L ' ' J 

m 


.3225 


.355 


PdV 


where 3 = 42.75 kb and ? = 8.390. Calculations have now been 
performed along the isotherm to 42.5 kb and 0.3225 gm/cc. 


The region from 42.5 kb to 80 kb along the 300°K isotherm is 

the region of complete transformation to coesite. Because 

little experimental evidence of the isothermal compression of 

coesite exists, we assumed a Murnaghan form using the same 

4 as quartz and replacing 3 with the isothermal bulk modulus 

c 

of quartz. The total energy is now: 
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The specific heat capacities at constant pressure for a-quartz 
and coesite in Region III are: 


C 

P 


a-quartz 


-.1324 + .0448 T 


3.074 X 10 

(200°K to 750°K) 


C 

P 


3-quartz 


14.65 + .00177 T 


(750°K to 1140°K) 


and 


C 

P 


coesite 


3.190 + .0297 T - 1.614 X 10 ^T^ 

(200°K to 920°K) 


The two equations for quartz are necessary to bridge the 
a-3 quartz transition at about 850°K. A linear mixing law 
between quartz and coesite is assumed in this region. 


C 

P 


coesite 


15.8115 + 1.132 X 10 ^T 


(920°K to 3000°K) 


The material transforms from 100% quartz at the beginning of 
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Region III to 100% coesite at the end of Region III. An aver- 
age specific heat capacity at constant pressure is calculated 
at each (3P/3V) interval of the calculation by assuinining appro- 
priate mixtures of the two C 's. AC is then calculated and 

P V 

used in the energy partition equation to calculate shock temper- 
ature . 

Regions IV and VA are the mixed regions for the transformation 
from coesite to stishovite. A linear mixing law for the corre- 
sponding specific heats at constant pressure results in tempera- 
tures which we consider too low. Thus Region IV was assumed to 
represent the region where 50% of the material is transformed 
to stishovite. Region VA represents the region where the 
remaining 50% of the material is transformed to stishovite. 

The appropriate specific heat at constant pressure is shown 
in Figure 12. 

Likewise, the energy from isothermal compression is the sum 
of the previous energies plus the flat region from V = .2975 
to V = .2275 and a Murnagham fit to the stishovite region. 



V = .2225 

.2975 

. 2375 

.2225 



+ 

-PdV + 

-PdV + 


V = .39 

.39 

' .2975 J 

.2375 


• .2225 

- “ 


/ \ 


''o/v 

. .2275 

> / 


We used the isothermal bulk modulus (3 ) and its first deri- 

s 

vative (5^) stishovite as given by Ahrens, et al. (19 70) 


34 



to fit the simplified Murnagham equation to the stishovite 
region rather than the Birch-Murnagham form with higher order 
corrections as given by Ahrens, et al. (1970). This was done 
to simplify the integrations and because any differences be- 
tween the two approaches seemed small compared with approxi- 
mations and estimations already made to reach this point in 
the calculations. 

Region VB is the region where material has been completely 
transformed to stishovite. The Cp and used in the temper- 
ature calculations is the same as previously given. 

Similar calculations are made using the solid quartz hugoniot. 
The various regions of the hugoniot and the 300°K isotherm are 
shown in Figure 13 and correspondingly in Figure 8. A discus- 
sion of the calculations will not be made but follows along the 
same lines as for the sandstone. Several comments are in order 
about the specific form of the solid quartz hugoniot. First, 
the HEL is taken as 45 kb which is close to an asmyptotic 
value of the HEL decay that is observed in solid quartz. This 
is interpreted by us to be a shock steady state equilibrium 
condition for the onset of the a-quartz/mixed region. Second, 
calculations are made for a-quartz transforming, directly to 
stishovite (B to C' in Figure 14) as well as the a-quartz-coe- 
site-stishovite (B to C) as was done for sandstone. The 
a-quartz to stishovite seems a more probable transformation 
within the shock front of laboratory experiments. The a-quartz- 
coesi te-s tishovite transitions are used because they appear to 
represent equilibrium conditions along the whole hugoniot. 

These two possible transformation processes should bound the 
actual transformation temperature conditions. Thirdly, because 
of the high resolved shear stress in a-quartz, a mean normal 
stress is computed through Poisson's ratio and the hugoniot 
stress . 
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Figure 14 High St 
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CONCLUSIONS 


Figure 14 illustrates the high strain rate phase boundaries 
for a-quartz-coesite and coesite-stishovite concluded from 
our interpretation of the hugoniot data and estimated iso- 
therms. The doubled lines represent the completed transition 
under shock loading and the dash lines represent the hydrostatic 
transitions. Error estimates are also shown. Lettered and 
numbered points in the figure correspond to discontinuities 
in the hugoniots shown in Figure 8. Discontinuities in the 
hugoniot temperature curves result from the simplifications 
of the hugoniots and isotherms . 

An estimate of the range of strain rates obtained in our shock 

2 11 

wave experiments is from 3 X 10 to 1.3 X 10 /sec. Strain 

rates for hydrostatic experiments are normally in the range 
— 6 ^3 

from 10 to 10 /sec. An obvious effect of the 14 to 17 orders 
of magnitude increase in strain rate on the polymorphism of 
quartz is a significant shift in the pressure-temperature plane 
of the phase boundaries of coesite and stishovite. 

Shock wave studies of quartz have application in at least two 
areas of geophysics. The first concerns the thermodynamic 
states of minerals present in the earth's mantel. Accurately 
measured P-V states assist in more accurate determinations of 
the thermodynamic properties and density states of quartz poly- 
morphs and thus help in determining the mineral assemblages that 
might be possible in the earth's mantel. 

The second area concerns meteorite craters and interpretations 
of the physical conditions that existed at impact. Many present 
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studies of such impacts base description of the P-V-T states 
on hydrostatic phase diagrams . The change in phase boundaries 
with strain rate as illustrated here might lead to a reevalua- 
tion of such interpretations. Because of the effect of strain 
rate as shown in Figure 14, an upper and lower limit may be 
established for the stability of quartz in the impact process. 
Previous interpretations of quartz should be viewed as only a 
lower bound. 
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APPENDIX 


HUGONIOT STANDARDS 


It has been found practical to induce ultra-high stresses 
by causing a high density metal plate to strike a target 
of the material of interest at high velocities. The basis 
of the present technique using a light gas gun as a projector 
is to impact sandstone targets with a flat plate of OFHC 
copper or Fansteel-77 (a tungsten alloy) . The hugoniot state 
in the target is measured by recording the impactor velocity 
prior to impact and the shock velocity induced in the speci- 
men. Computation of the resultant volume or particle velo- 
city requires, in addition to the above measurements, know- 
ledge of a hugoniot equation of state of the impactor material 
(McQueen, et al. , 1967). 

The hugoniot equations for the metals used in this program 
have been measured in this laboratory and are listed below. 

OFHC Copper : 

U = 3.964 + 1.463 u 
s p 

RMS Deviation of U = 

s 

Fansteel-77 : 

U = 4.008 + 1.262 u (km/sec) 
s P 

RMS Deviation of U = ± .021 km/sec. 

s 

A more complete description of the measurements on these 
metals may be found in Isbell, Shipman and Jones (1968) . 


(km/sec) 

± .009 km/sec 
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